SPATIOTEMPORAL CHAOS OF SOLITON IN A GENERALIZED 
SKYRME MODEL WITH THE MODIFIED SYMMETRY - 

BREAKING TERM . 



Nguyen Vien Tho [] 

Hue University, Hue, Vietnam 



Phu Chi HoaQ 



Dalat University, Dalat, Vietnam 



Abstract 



The chiral symmetry-breaking term of the Skyrme model with massive 
pion is modified to obtain the hedgehog profile function which is in best co- 
incidence with the kink-like profile function. For the modified Lagrangian, 
the minimum of the energy of the B=2 twisty skyrmion configuration is lower 
than the values for both the cases of the Skyrme Lagrangian with and with- 
out the non-modified symmetry-breaking term. The equations of motion for 
the time-dependent hedgehog of this model and for a generalizated Skyrme 
model including sixth-order stabilizing term are derived and integrated num- 
merically. The time evolution of soliton is obtained. We have observed the 
seft-exitation of soliton because of the fast developement of fluctuation. 

I. INTRODUCTION 

It is known that one can give a description of low energy hadron physics on 
the base of a semiclassical quatization of the soliton solution of the Skyrme model 
(skyrmion) [1 — 3] . The skyrmion configuration, denoted by f/(~af), is a map: R 3 — > 
SU(2), with the condition U — *■ 1 as |"a?| — > oo, required for finite energy . One 
could analyze in details the case of the spherically symmetry configuration, when 
the matrix field U(H?) has the form (the hedgehog): 



where r = and f is the Pauli matrices. Here F(r) is a profile function with 
the boundary conditions F(0) = nil (n is an integer) and F(r) — > as r^ oo. F(r) 
obeys a nonlinear differential equation which could be solve numerically. 

An analytic form of the skyrmion profile function seems to be very useful for 
many purposes. Attemps were made to obtain a such analytic approximation: (i) by 
computing the holonomy of SU (2) Yang- Mills instantons in R A along lines paralell 
to the time -axis [4]; (ii) by identifying the skyrmion profile function with the sin- 
Gordon kink field (if we replace x by r) [5]. The approach (ii) is very attractive 
because of its simplicity. Moreover, unlike the instanton approach, the sin-Gordon 
kink has a fixed scale, so there are no abitrary scale parameters which have to be fixed 
by hand in order to minimize the energy. With the explicit expression of the kink- 
like profile function F(r), in Sec. 2, the symmetry-breaking term of the Lagrangian 
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is modified to obtain the hedgehog profile function for B=l skyrmion which is in 
best coincidence with the kink-like profile function. Then we have calculated the 
energy of four different B=2 skyrmion configurations and lead to the result: the 
lowest value of the energy is obtained in the case of twisty configuration for the 
modified Lagrangian [6]. 

The results on numerical integration of the equation of motion for the time- 
dependent hedgehog indicate a dynamical chaostic character of fluctuations around 
the static soliton solution [7]. Such a study is interesting from the viewpoint of 
considering the Skyrme model as a nonlinear dynamical system. In Sec. 3, we find 
the equation of motion for the time-dependent hedgehog of the considered model. 
The results of numerical integration of the obtained equation are presented. The 
fluctuation of the profile function 8F(x, t) as well as the time dependence of the 
amplitudes of different modes of the fluctuations are plotted. 

The Skyrme's Lagrangian [1] is of fourth order in field derivatives. Various 
alternative models have been proposed which preserved the form of the original La- 
grangian while extending it to higher orders [8-12]. The incorporation of the higher 
order terms, on the one hand, improves the fit of observables, and, on the other 
hand, gives a reasonable physical interpretation for stabilizing terms in Lagrngian. 
For example, one could introduce a sixth-order term [8-9] 



£(6) = - £ 1 B »B„, (1.2) 



where is the baryon current 



p.fiua/3 

This term may be understood an representing effects of w mesons, while the Skyrme's 
fourth-order term should be viewed as representing effects of p mesons. In Sec. 4, 
we consider the dynamical behavior of the model which the Skyrme Lagrangian 
added by the sixth-order term (1.2 ) and the modified symmetry-breaking term. 
The discussion of the results is given in Conclusion. 

II. THE SIN-GORDON KINK FIELD AND THE CHOICE OF THE 
CHIRAL SYMMETRY BREAKING TERM 

The Skyrme's Lagrangian density takes the form [2,3] 

C = ^Tr (d»Ud»U + ) + ^Tr [{d,U)U + , (d u U)U + f . (2.1) 
One considers also the Skyrme's Lagrangian with a pion mass term [3] 

C' = C + C SB , (2.2) 

Csb = \m\FlTr (U + U+-2). (2.3) 



From (2.1), based on equations of motion, one found the nonlinear differential 
equation for the B=l hedgehog profile function 



xF' 
— + 



2 

^- + 2Sin 2 F 



,, „ f a2 Sin(2F) Sin 2 FSin(2F) . . 

F + Sin2F (F ) i—L -^-^ = 0, (2.4) 
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where x = eF w r is the dimensionless radial distance. 
The kink-like function [5] has the form 

F (x) = Aarctan (e~ x ) , (2.5) 

which satisfies the same boundary conditions. The kink-like profile function (2.5) 
has an exponential decay for large x [5]. The same of an asymptotic behaviour is also 
obtained when the Skyrme Lagrangian (2.1) is added by a chiral symmetry-breaking 
term as (2.3) [13]. Accordingly, the right-hand side of (2.4) is not zero, but equals 
to 

B 2 

'^-x 2 SmF, (2.6) 

where (3 = m^=140 MeV, e=4.84, F n =108 MeV. 

However, if one substitutes the kink-like profile function (2.5) in the left-hand 
side of (2.4) and compare with (2.6) (in Fig.la and Fig. lb), it is seen that they 
are different. In order to make the equality to be satisfied approximately, we must 
modify the pion mass term as following 
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3.5 x 10~ 7 — x 2 SinF. (2.7) 
4 

After modification the pion mass term as (2.7), we have plotted it in Fig.lc. We 
see that the approximate equality may be acceptable. So, one should be able to 
choice the symmetry-breaking term as follow [6] 

4T* = \m\FlTr (U + U+ - 2) , (2.8) 

where e = 3.5 x 10~ 7 . 

Now, we consider an ansazt has the form [13] 

N = {Cosk(f>Sin-&, Sink^Sinti, Costi}, (2.9) 

with this ansazt the soliton mass is 

M = M 2 + M 4 , (2.10) 

M 2 = 1 - r dxx 2 r dtfSm${(F') 2 + \k 2 + l}^4^}, (2.11) 
4 jo Jo x 2 

r°° , 9 r ,„„ „ rr Sm 2 l, 1l/w , 2 Sin 2 F~^Sin 2 F 

where 7 = lL ^r- l k is integer- k=l corresponds to the case of the sherically symmetry 
hedgehod, k> 2 to the case of twisty skyrmion configuration. The variation of (2.10) 
in F(x) give the following equation 

[x 2 + 2a sin 2 F]F" + 2xF' + [a(F') 2 - | - 26^^] sin(2F) = 0, (2.13) 



where 



a= r[A; 2 + l]sin», (2. 14) 

JO 
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b = k 2 f* smtidti, (2.15) 

JO 

which for k = 2 we have a = 10 and b = 8. 

From above-mentioned results and the formulas (2.11, 2.12), we consider the 
different cases of B=2 skyrmion configurations [6]. 

The first way of obtaining B=2 Skyrme field is to alter the boundary condition 
on the hedgehog profile function so that -F(O) = 2tt [14]. One can generate an 
approximation to this Skyrme field by using the kink-like profile function 

F (x) = 8arctan (V x ) , (2.16) 

substitute it to (2.11) and (2.12) with k=l, we get 

M a = M 2 + Mi~ 46.77. (2-17) 

The second way of obtaining B=2 skyrme field is to leave the boundary condition 
on the profile unchanged but to have the skyrme field rotate twice as rapidly as the 
radial vector x under a change in the azimuthal angle around an axis [15]. From(2.1), 
(2.11) and (2.12) with k=2, we get 

M b = M 2 + M 4 ~ 387. (2.18) 

In the third case, based on the Lagrangian with £$b (2.3) and k = 2,we obtained 
the equation of motion which is the equation (2.13) added by the term (2.6) in the 
right hand side. Solve this equation, we have nummerical dat file. Substituting it 
to (2.11) and (2.12) we get 

M c = M 2 + M 4 ~ 38.4 7 . (2.19) 

Analogously, in the last case, based on the Lagrangian with C^ d '^ and k = 2 
we obtained the equation of motion which is the equation (2.13) added by the term 
(2.7) in the right hand side. We find numerically F(x) (dat file) and substitute it to 
(2.11) and (2.12) we get 

M d = M 2 + M 4 ~ 27.57. (2.20) 



III. THE EQUATION FOR THE TIME -DEPENDENT HEDGEHOG 
AND THE TIME EVOLUTION OF THE SOLITON SOLUTION . 

We consider the Skyrme Lagrangian (2.1) is added by a chiral symmetry-breaking 
term (2.8). It is convenient to parametrize the SU(2) matrix field U(x) by the pion 
field iso vector if (x) [16,17] 

V(x)= \ + i Zf\Z , (3.1) 
l — ii] (x) r 



where 



t(*) = ^, (3-2) 
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and t — (ri,r 2 ,T3) are the Pauli matrices. In the parametrization (3.1) the expres- 
sions of the Cartan forms (d^U)U + for the SU(2) group have been calculated [16,17], 
and on the base of these expressions we get the expressions of the Lagrangian (2.1) 
an (2.8) in the form 



C = 



F 2 (d^Jfd^Jf) 4 [dfj/rf x d^lf]' 



(l + lf 



(l + lf 



f(mod.' 
'-SB 



~mlF 2 



2 * ^(l + t 2 )' 
The time- dependent hedgehog corresponds to the following ansatz 

'F(r,t) 



if (r, t) = tan 



r 
r 



(3.3) 
(3.4) 



(3.5) 



where F(r,t) is the profile function. Corresponding to this ansatz, C and C^ d ^ are 
given by 



F 2 ,-. , 2sin 2 F, 1 sin 2 F , , , sin 2 

£ = -^i(F) 2 -(F) 2 - + ^{2(F) 2 - 2(F ) 2 - — — }, (3.6) 



r 2 2e 2 r 2 r 



The Hamiltonian is given by 



a ( 2 , r^ 2 r/^2 /n 'x 2 2 sin 2 F , 1 sin 2 F_.-. 2 ,_,/,. , sin 2 

H = 4n r 2 dr{-^{(F) 2 + (F) 2 + — } + ^^{2(F) 2 + 2 F ) 2 + ~^}+ 

J 8 r z 2e z r z r z 

£m 2 F 2 
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The variational equation for the profile function is 



1-cosF)}. (3.8) 



xF' 



+ 



x 



+ 2Sin 2 F 



>? - F) + Sm2F(f — F 2 ) — S -^l - ^ FSm < 2F > 



2T 



e(3 2 



-x 



2 SinF = 0, 



(3.9) 



where (3 = zjfc, m^=140 MeV, e=4.84, F w =108 MeV [ 14 ] , and x = eF n r, 
r = eF w t are the dimensionless distance and time. The primes and the dots mean 
the derivatives with respect to x and r, respectively. Hereafter t is always understood 
as the dimensionless time. 

It is convenient to write F(x,t) in the form 



F(x,t) = F(x) + 5F (x,t) , 



(3.10) 
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where F(x) is the profile function of the static hedgehog. To find 5F (x, t) we use 
the following boundary conditions 



5F (0,t) = 5F(L,t) = 0. 
This condition is automatically satisfied by the harmonic expansion 



(3.11) 



N-l 



JTIX 



) 



5F(x,t) = Aj(t)wa. 



L 



(3.12) 



where L is the size of the spatial volume, N is the number of the points of the 
discretized spatial variable x, Aj(t) are the amplitudes of jth fluctuation modes. 
These amplitudes are obtained by inverting the series given in (3.12) 



In our calculation we choose L=16, N=128 and the initial excitation mode is j 
= 16. We denote A = Ai 6 (0). We have studied the time evolution of the system 
for the pertubation parameter A = 0.1. The equation (3.9) should be reduced to 
a system coupled second order differential equations for the time variable t. We 
solve this system by using the Runge-Kutta procedure and obtain the solutions for 
5F(x,t) and A 6 (t). We have plotted the fluctuation 8F (x,t) at t=0, 100, 200, 300 
and 500 in Fig. 2a to Fig.2e, respectively. We see that apart from large fluctuations 
near x — 0, \8F(x,t)\ ~ 0.1 . This is understandable as the Skyrmion dynamics 
is dominated by the small x-region and the deviation from the Skyrmion is small. 
The dynamical behavior of the system is understood better by observing the time 
evolution of the amplitudes of various harmonic modes. In Fig. 3a to Fig. 3c we have 
plotted Aj(t) for j=8, 16, 32 , respectively. It is seen that the amplitude of the initial 
mode decreases gradually and it has some kind of periodicity in the variation, while 
amplitudes of other modes increase on the average. The time evolution of the mode 
of j=8 and j =16 has a periodic behavior when it was consider in a small interval 
of the time. It is to be expected that after a longer interval of time, all the modes 
would be of the same magnitude leading to "thermalization" and "spatio-temporal 
chaos" . 

IV. DYNAMICAL BEHAVIOR OF SKYRMION IN A GENERALIZED 
SKYRME MODEL . 

We consider the Lagrangian 



where + is the Lagrangian of the Skyrme model (2.1), is given by 
(1.2 ), £^ od '- ) is given by (2.8). In the parametrization (3.1), the expressions of the 
Cartan forms (d l JJ)U + for the SU(2) group have been calculated [16,17]. We get 
the Lagrangian (1.2) in the form 




(3.13) 



(mod.) 



(4.1) 



SB 




<-{(dJf) 6 - 3(dJf) 2 (d a JfdPjf)(dpJfd a Jf) + 



6 



2(d u Jfdi 3 Tf)(d f3 Jfd a Jf)(d a Jfd u Jf)}. (4.2) 

By the same caculation in Sec. Ill, we have obtained the variational equation for 
the profile function 



, x 



2 „, QA„i 



C _ + 2S ,„>F + \^)(F - F') + (Sin (2F) + I ^W) ^ _ F *,_ 

2 4 ar 4 x 2 4 

where x = eF n r, r = eF n t are the dimensionless distance and time. The primes 

F 2 e 2 e i 

and the dots mean the derivatives with respect to x and r, respectively, 7 = * | , 
/3 = We choose £g = 5/m 2 , and t is understood as the dimensionless time. 

Based on the expressions (3.10) to (3.13) and the choice of L=16, N=128, A = 
-^16(0) = 0.1, we consider the dynamical behavior of the system with Lagragian ( 
4.1) by observing the time evolution of the amplitudes of various harmonic modes. 
We have plotted the profile function of the static hedgehog and the fluctuation 
5F (x,t) at t = 0, 100, 200, 300 in Fig. 4a to Fig. 4d. It is a clear indication that 
SF (x, t) has violent fluctuations for substantially longer intervalls in x. 

The plots of Aj(t) for the cases of j= 8, 16, 64, 127 ( Fig. 5a to Fig. 5d ) indicate 
that the fluctuation amplitudes in the generalized Skyrme model which Lagrangian 
added by the sixth-order term develop much faster than in the Skyrme model. For 
this model, the amplitudes of "spontaneous" fluctuations appear after t as 150. One 
can say that a self-exit at ion of soliton takes place after t « 150. 

V. CONCLUSION 

We have calculated the energy of four B=2 skyrmion configurations. The com- 
parison of the values of the energy in four cases shows that in the case of modified 
Lagrangian one could obtain B=2 configuration with lowest energy. That is, this 
configuration is more close to the real energy minimum of B=2 skyrmion. 

The B=l hedgehog kink-like profile function (2.5) is a convenient analytic ap- 
proximation to the numerical solution. But how can modify the Lagrangian to 
obtain the profile function (2.5), the modification made in this paper is one of the 
answers to the question. 

By integrating the equation (3.1) for the time-dependent hedgehog we have ob- 
tained the information about the dynamical behavior of the sotiton in the model 
(3.3, 3.4). We have plotted 5F(x,t) at various moments and the development of 
amplitudes of fluctuation modes Aj(t). The plots of Aj(t) for the case A=0.1 in- 
dicate that the process of thermalization takes place sooner, and the fluctuation 
amplitudes in the considered model develop much faster than in the Skyrme model 
(see [7]). Besides, in [7], for A=0.1 the amplitude of the initial (32nd) mode has 
extremely regular periodic behavior while in our initial (16nd) mode the amplitudes 
decreases gradually on the average. 

When higher order terms are included, the theory becomes much more nonlinear, 
so it is clearly that the phase space around the soliton solution in the generalized 
Skyrme model was shown to be more stochastic than in the original Skyrme model. 
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